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Abstract 

The X-ray active RS CVn binary system a Gemino- 
rum was observed during a long-duration flare with XMM- 
Newton. We compare line emission during the flare with 
line emission measured from a previous Chandra observa- 
tion during quiescence, in a model independent way. We 
flnd that in addition to an overall 25% flux increase, which 
can be ascribed to variations in the system's quiescence ac- 
tivity over the 15 months that passed between the obser- 
vations, there is a hot plasma component of kTe > 1 keV 
that arises with the flare. The hot component is mani- 
fested primarily by high charge states of Fe and by a vast 
continuum, but with no temperature or abundance effect 
on the cooler {kT^ < 1 keV) plasma, except perhaps for 
Mg that is anomalously enhanced. 
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1. Introduction 

The interplay between steady coronal emission and coro- 
nal flares has been a subject of ongoing research for many 
years. X-ray line resolved spectra available with Chandra 
and XMM-Newton allows now for unprecedented plasma 
diagnostics, one of which is elemental abundance measure- 
ments. 

The RS CVn binary system a Gem (HD 62044, HR 
2973, HIP 37 629) is X -ray bright (logLx ~ 31.0 ± 0.2 
erg s~^; lYi et al. 1997|l and well observed in all wave- 
lengths. It has a rather long period of 19.6045 days (Duemm- 
ler et al. 1997) for such an X-ray bright star. The primary 
star is a Kl III type red giant. Little is known of the sec- 
ondary as it has not been detected directly at any wave- 
length, but restrictions to its mass and luminosity suggest 
that it is most likely a late type main-sequence star of 
one solar mass or less l|Duemmler et al. 1997|l . While most 
previous observations found it to be a relatively steady 
emitter, a very large flare has been detected in Decem- 
ber 1998 with EUVE l|Sanz-Forcada et al. 200^ . Another 
flare was detected in April 2001, both in X-ray and in radio 
using XMM-Newton and the VLA HGiidel et al. 2002|l . 



In this poster we present a direct comparison of the 
soft X-ray emission from the April 2001 flare with a qui- 
escent spectrum obtained in December 1999 in a model 
independent way. Our goal is to use only the line emis- 
sion for the comparison as the lines constrain the plasma 
properties more tightly than the continuum. 

2. Observations and data analysis 

2.1. The observations and light curves 

The target a Gem was observed by XMM-Newton with the 
Reflection Grating Spectrometers (RGS) on April 2001, 
for a total exposure time of 54 ks, of which we obtained 
data in the range of 6 - 38 A. Chandra observed the target 
on December 1999 for 100 ks with the Low Energy Trans- 
mission Grating (LETG) spectrometer in the Continuous 
Clocking (CC) mode for which reliable data were available 
only for the 2 - 20 A range. 

In the long duration flare observed with XMM-Newton, 
a Neupert effect was detected using the EPIC-PN camera 
on board XMM-Newton and the VLA. The Neupert effect 
supports the scenario in which the radio emission is due 
to fast particles spiraling down the magnetic loop, emit- 
ting gyrosyncrotron radiation, and the X-ray emission is 
due to chromospheric plasma heated by these particles. 
Thus the radio flux depends on the instantaneous number 
of particles spiraling down while the X-ray flux depends 
on the total energy transferred to the plasma, and the 
approximate relation is -jf.Lx ck L Radio- These results, 
reported by Giidel et al. (2002), are shown in Figure ^ 
The RGS light curve 1st and 2nd order for both instru- 
ments, rebinned to 500 s bins, is presented in Figure |21 
The flare light curve can be contrasted with the flat light 
curve from the Chandra observation (using all orders of 
LETG), which is shown below in the same figure. Since 
the Chandra / LETG light curve is very flat we used this 
observation as a quiescence state reference. 

2.2. Comparison of flaring and quiescent spectra 

The fluxed spectra extracted from the instruments are 
presented in Figure |3 The overlapping waveband of the 
instruments is only 6 - 20 A. The most outstanding dif- 
ference is the large addition of continuum emission in the 
flare spectrum. The total flux in the overlapping region is 
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Figure 1. Light curves of a Gem,, oht,ained on April 2001, 
showing the Neupert effect, reproduced from Giidel et al. 
(2002). (a) XMM-Newton EPIC PN light curve, binned 
to 500 s. The sm,ooth solid curve has been obtained from 
a Chebychev polynomial fit of order 9. (b) Time deriva- 
tive of the smoothed X-ray light curve. Solid curve: Time 
derivative of the Chebychev polynomial fii above. Dotted 
curve: Derived from a boxcar- smoothed light curve (box- 
car length = 11; not shown). For the latter derivative, 
the larger error bar in the upper right corner illustrates 
the absolute uncertainty for any single derivative, whereas 
the smaller error bar indicates the relative scatter between 
nearest neighbors and is smaller due to correlations in- 
troduced by the smoothing, (c) VLA 6 cm light curve, 
binned to 750 s (observing scan length). The inset shows 
the cross-correlation function of the radio light curve and 
the X-ray time derivative, computed for the time interval 
marked with a horizontal bar above the second radio flare. 



increased by a factor of 2.35 in the flare spectrum, relative 
to the quiescent one. 

In order to compare line emission, we have taken into 
account the different Line Spread Functions (LSF) of the 
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Figure 2. Light curves for a Gem in time bins of 500 s. (a) 

April 2001 RGS (1 and 2 combined) observation (1st and 
2nd orders), (b) December 1999 LETG/ACIS observation 
(all orders). 
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Figure 3. a Gem spectra, (a) Flaring XMM-RGS spectrum, 
April 2001. Averaged spectrum of the two RGS instru- 
ments in 1st order, (b) Quiescent state Chandra-LETG 
spectrum, December 1999. Averaged spectrum of the two 
1st orders. 



instruments. LETG has approximately a Gaussian LSF 
with FWHM « 0.05 A, while RGS has a slightly wider 
LSF of FWHM K. 0.07 A and a more complex shape with 
a narrow peak. For the comparison we fold the LETG 
data through the RGS response, but in order to reproduce 
the RGS profile correctly, we first need to represent the 
LETG bright lines as delta functions with their correct 
measured flux. This is done by Gaussian fitting to the 
lines with a local continuum. We then subtract the fitted 
Gaussians from the spectrum to get a residual continuum. 
Finally we fold the delta-function lines plus the residual 
continuum (thus ensuring the correct flux is retained with 
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Figure 4- The LETG source spectrum compared with the 
fitted Gaussians + residual continuum spectrum. The 
delta-functions and continuum were passed hack through 
the LETG response and subtracted from the LETG spec- 
trum to give the plotted residual. Random fluctuations are 
smoothened by to the second convolution with the LSF. 
The residual shows a very good match, to the order of the 
random fluctuations, indicating the delta-functions -h con- 
tinuum and the source spectrum cannot be distinguished by 
the LETG spectrometer (and also by RGS). 

high precision even if there were some inaccuracies with 
the fitting) via the RGS response. 

As a consistency check, we also fold the processed spec- 
trum (delta-functions -t- residual continuum) back through 
the LETG response and compare to the original data. The 
results are presented in Figure^ and they show that the 
source spectrum and the processed one cannot be distin- 
guished by the LETG instrument. Note that in the pro- 
cessed spectrum, after being passed through the LETG 
response, the statistical fluctuations are smoothened out 
by the second convolution with the LSF. 

In order to investigate the change in line intensities 
from quiescence to flare we chose to scale the emission lines 
by the intensity ratio of the bright, well isolated, 0^+ Lya 
line at 19 A. Accordingly, the quiescent spectrum was mul- 
tiplied by a factor of 1.25. To illustrate the comparison, 
the remaining continuum difference was also added to the 
quiescent spectrum to match the flaring continuum. The 
end result of this process is shown in Figure where the 
flare and scaled quiescent spectra are plotted one on top 
of the other. Note that the continuum is the flare contin- 
uum for both spectra. Only the lines should be compared. 
This enables a comparison of the relative line intensities 
directly between the two phases. 

3. Results 

As seen in Figure |S1 the relative line intensities at wave- 
lengths longer than 12 A match very well in the two spec- 



tra. The fluctuations in the residual around the brightest 
lines at 19, at 13.5, and at 12 A are caused by slight inaccu- 
racies in the LSF of the RGS intensified by the large steep 
peak. The residuals however average to nearly zero and in 
the 19 A line, exactly zero, by design of the scaling. The 
lines above 12 A come from ions with maximum formation 
temperature of under kTf, 600 eV and maximum emis- 
sivity temperature (the temperature for which the emitted 
line is brightest) of under 650 eV. Lines at shorter wave- 
lengths of Fe^^+ and 81^^+ at 12.84 and 6.65 A, respec- 
tively, whose maximum emissivity temperature is 850 eV 
also match well in the two spectra. 

We conclude that the ions with maximum emissivity 
temperatures comparable to that of Fe^^^ 1 keV) and 
below were not affected by the flare, except their emis- 
sion measure increased by a constant factor of 1.25. This 
means that the plasma at temperatures below ~1 keV has 
the same abundance composition in the two observations 
and no change in the shape of the emission measure distri- 
bution is observed. However, the total emission measure 
did increase by the 1.25 factor by which we multiplied the 
quiescent lines. This means that either the flare is evapo- 
rating plasma from the chromosphere with the same chem- 
ical and thermal composition as the quiescent corona or, 
more likely, the 1.25 enhancement is unrelated to the flare 
and the amount of relatively cold corona increased by this 
factor over the 15 months that passed between the obser- 
vations. Previous observations with ROSAT spanning 1.5 
yrs. showed changes of the same order l|Yi et al. 1997|l . 
This might be related to the system's activity cycle or 
changes in active regions and not to the flare. 

More evident changes in the lines are observed for 
highly ionized iron lines of Fc^"+ to Fe^'^^ below 12 A, 
which correspond to temperatures of 1 keV and above. 
Since the relative intensities of the lower-T Fe lines do 
not change during the flare, this is clearly a temperature 
effect rather than an abundance effect. This and the en- 
hanced continuum at short wavelengths without change 
to the lines, indicate the appearance of a very hot com- 
ponent with little effect on the cooler plasma. The only 
other ions that probe this temperature range and appear 
in the RGS spectrum are Al^^+ at 7.17 A, which is too 
weak to be used and 81"'^^+ at 6.18 A, which does seem to 
be enhanced, but is just at the edge of the RGS spectral 
range and, thus, is not clear enough to be used. 

Mg-'^^+ is the exception to the rule. Although its max- 
imum emissivity temperature is 840 eV, it does show a 
large increase in the Lya line at 8.42 A line, by a fac- 
tor of 2.3 ± 0.5 beyond the overall 25% enhancement. 
On the other hand, the other ions with similar maximum 
emissivity temperatures, namely Fe-^^+ and 8i-^^+ (850 eV 
and 860 eV) do not show the same effect. Since Mg-^^+ 
is a Hydrogen-like ion, its emissivity drops off slowly be- 
yond the 840 eV peak and extends to higher temperatures 
than those of Fe^^+ and Si^^"*". Thus, we could be see- 
ing its emission originating at higher temperatures in the 
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Figure 5. Line intesities of the flaring and quiescent states compared. The quiescent lines are multiplied by a factor of 
1.25 and an ad-hoc continuum is added to match the flaring continuum. The plotted residual is the quiescent curve 
subtracted from the flaring spectrum, with corresponding error bars. 



flare. The fact that no change is detected in the Mg^*'+ 
lines which are emitted primarily at lower temperatures 
of ~540 eV reinforces our conclusion that the cold gas is 
not affected by the flare. However, since we concede that 
ascribing Mg^'^^ to the flare is somewhat ambiguous, an 
abundance enhancement of Mg cannot be totally ruled 
out. 



4. Conclusions 

1. The 25% increase in flux between 1999 and the flare 
of 2001 may be due to gradual variations in the soft 
X-ray activity rather than a result of the 2001 flare. 

2. The flare shows a two component behaviour: cooler 
plasma which is not effected by the flare in a detectable 
amount and plasma at kT^ > 1 keV which is signifi- 
cantly enhanced. The flare is manifested unambigu- 
ously by a strong continuum that rises down to the 
RGS limit at 6A (while the quiescent continuum peaks 
around ~ 13 A) indicating temperatures of kT^ ^ 
2 keV, higher than the temperatures producing the 
lines. This suggests that the large flare is localized and 
does not disturb the surrounding cooler plasma. 

3. No low-T evaporation nor enrichment of the cooler 
parts of the corona is detected. The flare did not change 
the relative abundance of elements in the corona in a 
detectable amount, with the exception of maybe Mg. 
Additional ions of Si, S, Ar, and Ca that could indicate 



changes at very high temperatures fall outside of the 
RGS band (i.e., under 6 A). In relation to the chro- 
mospheric evaporation scenario, this suggests that the 
amount of evaporated material is small relative to the 
quiescent corona. As the hot evaporated plasma rises 
from the denser chromosphere to the higher corona it 
expands and its emission measure drops proportionally 
to the electron density. Hence, the effect on the emis- 
sion measure of the cooler coronal plasma is negligible 
and remains undetected. 
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